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Abstract 

A new type of (3 model of intracluster gas distribution is established, which is well fitted 
to the observed X-ray surface brightness profile except for outermost region. It employs the 
exact King model of the self-gravitating system of collisionless particles instead of the widely 
used approximated analytical formula for the dark matter distribution. The small difference 
between the modelled and the observed X-ray surface brightness profiles would disappear if the 
infall of matter from the intercluster space is taken into account. The number of parameters 
of our model is the same as for the isothermal f3 model, but the agreement between our model 
and data points is satisfactory even when the fitting with the isothermal (3 model remains a 
large central excess. In addition, the enhanced central mass concentration in our model is high 
enough to resolve the discrepancy between the X-ray and the lensing measurements of cluster 
mass. These facts suggest that most of the central excess might be no more than a ghost em- 
bossed by an inadequate fitting. Therefore the central excess would not necessarily be linked to 
the existence of cooling flow. It is also shown that the obtained values of f3 (our model) ~ 0.5 
support the theory of non-gravitational intracluster gas heatings comparable to those reached 
through gravitational collapse. On the other hand, the value of (3 ~ 0.6 obtained for the 
isothermal (3 model is explained as showing the scarcity of matter in the space surrounding 
each cluster, brobably due to the infall of matter towards it over a long time. 

Key words: X-rays, cluster of galaxies, intracluster gas temperature, density distribution, 
dark matter 



1. Introduction 



The standard model of mass distribution in clusters of galaxies assumes that both gas and 
dark matter are in hydrostatic equilibrium within the binding cluster potential. Recently many 
studies have compared mass measurements for clusters of galaxies using the X-ray and the 
lensing techniques. Miralda-Esqude & Babul (1995) reported a discrepancy by a factor of ~ 2 
between the X-ray and the strong-lensing mass determinations for A1689 and A2218. Wu & 
Fang (1997) concluded that the mass measurements from the velocity dispersion of optical 
galaxies and the strong-lensing data were generally in good agreement but were typically a 
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factor 2-3 larger than those from the X-ray data. In contrast to this, the comparison of the 
weak lensing and the X-ray mass measurements at cluster's outer regions has inferred good 
agreement between them (e.g. Squires et al. 1996; Smail et al. 1997). These facts show that 
the isothermal (3 model widely used in the X-ray data analyses is useful preferably only in the 
outer region of each cluster. 

It is therefore inebitable to modify the model of intracluster gas distribution in order to 
resolve the discrepancy between the X-ray and the strong-lensing mass measurements. The 
double f3 model of Xue & Wu (2000) is a possible one. However, as shown by themselves, the 
correction by the double (3 model is still insufficient to resolve the discrepancy. In addition, 
the physical foundation of the double (3 model is unclear, and the mathematical formulas to 
reconstruct the actual form of gas distribution are fairly tangled and unphysical. For example, 
the restored gas distribution is not a simple superposition of two separated components. Models 
of intracluster gas distribution based on the NFW model of dark matter halo having central cusp 
(Navarro, Frenk & White 1996, 1997; Moore et al. 1999) or on the purely phenomenological 
dark matter density distribution with no central cusp (Burkert 1995) are suitable only for 
individually limited types of clusters. Such a state of things demands us to establish a new 
powerfull model of intracluster gas distribution, on the basis of firm physical foundations, which 
can resolve these problems all at once. 

In this paper, we propose a new type of (3 model of intracluster gas distribution in Section 
2. Then we determine the parameter values of our model best fitted to observed X-ray surface 
brightness profiles of three clusters of galaxies and compare them with the results obtained by 
using the conventional isothermal /3 model in Section 3. The good agreement between the X-ray 
and the gravitational lensing mass measurements based on our new model is also shown there. 
In Section 4 we discuss the new features of clusters of galaxies on the basis of our model. Brief 
comments are also given on the works based on other influential models. Finally, the present 
results are summarized in Section 5. 

Instead of the recent WMAP results showing the existence of non-zero A = 0.73 ± 0.04) 
and the value of H = 72 ± 5kms _1 Mpc _1 (e.g. Bennett et al. 2003; Spergel et al. 2003), we 
assume an Einstein-de Sitter (Qm = 1.0, Q\ = 0.0) model of the universe and adopt a Hubble 
constant of H = 50kms _1 Mpc _1 in this paper, except for Section 4, in order to make it easy 
to compare the present results with many previous works which used those old cosmological 
parameter values. 

2. Theoretical Foundations 

Here we outline the theoretical basis for making up the model function of the intracluster gas 
distribution and check in detail the general properties of the conventional isothermal (3 model. 
Then we establish a new type of model of intracluster gas distribution, which is expected to be 
more physical and more practical than the isothermal (3 model. 

2.1 Basic relations 

The standard model of mass distribution in clusters of galaxies assumes that both gas and 
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dark matter are in hydrostatic equilibrium within the cluster potential. That is, 

1 dP(r) _ G MW ; (1) 



p g (r) dr r 

where G is the gravitational constant, r is the three-dimensional distance from the cluster 
center, and p g (r), P(r) and M(r) are the intracluster gas density, pressure, and the total 
mass (including dark matter) within the radius r, respectively. The equation of state for the 
intracluster gas is 

where mn, A; B , p and T g (r) are the hydrogen mass, the Boltzmann constant, the mean molecular 
weight and the temperature of intracluster gas at r, respectively. From Equations (1) and (2) 
we obtain 

/c B T g (r)r f dha.p g (r) i d\nT g (r) 



G/jm H \ amr alnr J 

This expression is available in good approximation even when a radial flow of intracluster 
gas exists, so long as the flow velocity is much less than the thermal velocity of gas particles 
(Miyoshi 1994). 

On the other hand, the collisionless dark matter particles satisfy the following relation (e.g. 
Binny & Tremaine 1994, hereafter BT) 

1 <WW(r)) + 2 / _ I N = _ G MM_ (4) 



Pdm(^) dr r \ r 2 / r 

where of (r) and of (r) are the radial and the tangential (perpendicular to the radial component) 
velocity dispersions, and Pdm(^) is the dark matter density. Equation (4) gives 

M(r) = -^P ( + + 2A(r)) , (5) 

G \ alnr amr / 

where A(r) = 1 — of(r)/2af (r) indicates the degree of anisotropy in the velocity distribution 
of dark matter. From Equations (3) and (5) we can introduce a parameter (3 defined by 

2 d\np g riTnT g 

ft — <Jr — dlnr dlnr /g\ 

pm-R dlnr dlnr 

When the intracluster gas is isothermal and the dark matter velocity dispersion is isotropic and 
r-independent, that is, 

f/ln7; c/hi a,. ^ _ q ^ 



Equation (6) is reduced to 



dlnr dlnr 



dlnr dlnr 



This is equivalent to 



(9) 



It should be noted here that Expression (9) is still valid in good approximation even if A, T g 
and ay are r-dependent, so long as the following inequalities hold 



dlnTrr 



dlnr 
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(10) 



2.2 Isothermal (3 model 

When the density distribution of dark matter is assumed to have the same form as the 
asymptotic density distribution of self-gravitaing isothermal gas sphere, 



p DM (r)oc[l + (r/a) 2 r 3 / 2 , 



(11) 



which is fairly well fitted to the galaxy distribution in clusters of galaxies (King 1972), Expres- 
sion (9) gives the isothermal f3 model of intracluster gas distribution (Cavaliere & Fusco-Femiano 
1978; Jones & Forman 1984) 



p g (r)=p g (0)[l + (r/ay 



(12) 



where p g (0) is the central gas density and a is the core radius. If the intracluster gas is 
isothermal and the radial gas density distribution is of the form of Equation (12), the X-ray 
surface brightness distribution which is proportional to the integral of p g (r) along the line of 
sight has the functional form of 



S(R) = S(0)[1 + {R/af 



-3/3 



(13) 



where 5(0) is the central surface brightness and R is the projected radial distance from the 
cluster center on the sky plane. Many clusters, however, are not adequately described their X- 
ray surface brightness profiles by this model but show some central excesses. Fabian & Nulsen 
(1977) and Jones & Forman (1984) interpreted these central excesses as being due to the cooling 
flow of intracluster gas. The principle of their idea is as follows. Since the X-ray emissivity 
is proportional to y^p 2 , the lower gas temperature yields the higher X-ray emissivity, when 
the gas pressure T g p g remains unchanged. The last condition is required to keep the balance 
between the pressure gradient and the gravitational force. 

Here we have to remember that the density profile of self-gravitating isothermal gas sphere 
with a finite value at the center has a central hump in general (see e.g. BT). The X-ray 
surface brightness profile simulated using this solution also has a central hump, in contrast to 
the isothermal (3 model. It is therefore probable that the central excess in the X-ray surface 
brightness profile remained in the fitting with the isothermal (3 model might be mostly no 
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more than a ghost. Accordingly the central excess does not necessarily mean the existence 
of cooling flow. Besides, the cooling flow model cannot resolve the discrepancy between the 
lensing and the X-ray mass measurements mentioned in Section 1. Because, the cooling flow 
model requires only the enhancement of gas in the central region, and the total mass density 
distribution (mostly due to dark matter) remains almost unchanged. 

These facts strongly provoke us to make up a new model function of PdmM, and therefore 
p g (r) through Expression (9), which is expected to make clear the vagueness concerned with 
the central excess as well as to resolve the discrepancy between the X-ray and the lensing mass 
measurements. 



2.3 The new model we propose here 



Here we propose a new model function of p g (r) obtained from Expression (9) by using exact 
King models for Pdm(^), which is well fitted to the observed X-ray data as shown below. 

It is well known that the stellar distributions in globular clusters and elliptical galaxies are 
very well described by King models (Michie 1963; Michie & Bodenheimer 1963; King 1966, 
1981) defined by 



p(r 
P 



^ = e^W ( y*W/,)-^g)(l + ?|W), (14) 

where pi is a normalization factor, a is a constant of the order of the dispersion velocity, 
erf (x) is the error function defined by 

erf (x) = [ e~ t2 dt, (15) 
V 71 " Jo 

and \l/(r) is a solution of the following Poisson equation 



In terms of p = pj p , ^ = ^/o" 2 , and r = r/r , where p = p(0) and r = ^9a 2 /AnGpo 
(the King radius; see BT for details), Equations (14) and (16) are reduced to 

p(f)=X(*(f))/X(* ), (17) 

and 

r 2 dr y dr J X(tf„) 
where = ^(0), and X(x) is defined by 

X(x) = e x erf (Vx) - y/4x/ir(3 + 2x)/3. (19) 

The solution of Equation (18) gives a normalized density profile p(r). The boundary condition 
adopted here is d^/df = at f = 0. The density profiles thus obtained for finite values of 
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= 12, 9, 6, and 3 and their projected density profiles are shown in Fig. 4-9 of BT, each of 
which shows a central hump more noticeable for larger values of \^o- The dispersion velocities 
for these solutions are nearly constant in the central region and gradually decrease with r in 
the outer region (see Fig. 4-11 of BT). Expression (10), the condition required for the validity 
of Equation (9), is therefore well satisfied. 

We can construct a new model of intracluster gas density p g (r) from Equation (9) by sub- 
stituting the solution p(r) of Equation (14) for Pdm(?")- This model is necessarily expected to 
have a central hump and therefore to yield some central excess in the X-ray surface brightness 
distribution when fitted by the isothermal (3 model. The X-ray surface brightness is propor- 
tional to S{R), the integral of p 2 g {r) along the line of sight. In Fig. 1 we plot S{R) simulated for 
our models of \& = 12, 9, 6, and 3 in case of f3 — 1.0. The expected central humps are clearly 
seen for large values of In the limit of \j/ — ► °o, S(R) converges to one for the isothermal 
sphere, i.e. S(R) oc R~ 3 at large R. 

As was pointed out by Cavaliere & Fusco-Femiano (1978), the exact equilibrium solution 
can be obtained by solving the full Poisson equation containing the total (dark matter plus 
baryonic matter) density. In our model, however, the Poisson equation (16) or (18) includes 
dark matter only. That is, the gas (baryonic matter) density is not included in the right hand 
side of these equations. Gas particles are distributed in the gravitational potential made up 
by dark matter. In addition, there is no foundation on which the validity of Equation (7) is 
supported, and only the observational data suggest the validity of the conditions described 
in Expression (10). Therefore, Equation (9) should be still considered to be an approximated 
model of real intracluster gas distribution. Nevertheless, the agreement between our new model 
and the actual intracluster gas density profile is expected to be considerably improved in com- 
parison with the isothermal (3 model. Because, the fractional mass content of dark matter is 
overwhelming there and the gravitational field within each cluster is practically governed by 
dark matter. 

3. Comparison with Observations 

We present here the results of fitting our model to the observed X-ray data of three clusters 
of galaxies, A 383 (Smith et al. 2001), A 2163 (Elbaz, Arnaud & Bohringer 1995) and A 2390 
(Bohringer et al. 1998), compared with the results of fitting by the isothermal (5 model. A 383 
and A 2390 are classified in cD clusters, but A 2163 is not. This cluster is classified in the Abell 
richness class 2. Since giant arcs are observed in the central regions of these clusters, we can also 
compare their projected masses, M pro j, the integral of the mass density obtained from X-ray 
data along the line of sight within the arc radius, with their (strong) lensing masses, Mi ens . The 
strong-lensing masses of A 383 and A 2390 are obtained from Smith et al. (2001) and Pierre 
et al. (1966), respectively. For A 2163, the strong-lensing mass determined by assuming the 
arcs lie at the Einstein radius is 4 x 1O 13 M , but this value exceeds the weak-lensing masses 
obtained at some outer radii (see Fig. 2 of Squires et al. 1997). Therefore we adopt for Mi ens 
the value of weak-lensing mass estimated by extrapolation from outer radii. For the intra- 
cluster gas temperature of A 383, Smith et al. (2001) adopted the value of T g = 7.1 (±2) keV 
estimated from the cooling-flow-corrected cluster Tx-L Bo i relation obtained by Allen & Fabian 
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(1998). The use of thus obtained value of T g is, however, not suitable for our analysis because 
of the large error bars accompanied with the relation. The archived Chandra data of this clus- 
ter yield the emission-weighted mean temperature to be 3.5 ± 0.2 keV and 5.2+Qg keV for the 
central (Or < 15", where Or is the angular distance from the cluster centre) and surrounding 
(15" <0 R < 120") regions, respectively (Furuzawa 2002). We adopt the value of T g = 4.0 keV 
for the present analysis. The observed redshifts (z), the gas temperature (T g ), the projected 
distance of giant arc from the cluster center (i? arc ), and the lensing mass (Mi ens ) of each sample 
cluster, we adopted for the present analysis, are given in Table 1. 

In Figs. 2-4 we plot the results of fitting for our model and the isothermal (3 model to the 
observed X-ray surface brightness profiles of these clusters. The best-fit values of model pa- 
rameters and the resultant M pvo -J M\ ens ratios for the two models are summarized in Table 2, 
in which the projected mass obtained from the isothermal f3 model is indicated by MrT iPr oj- 
The agreement between our model with best-fit parameter values and the observations is very 
satisfactory for all the sample clusters aside from their outermost regions. The differences 
between the model curve and data points are much less than the statistical errors for A 383 
and A 2390. Thus our model is in general much more suitable to reproduce the X-ray surface 
brightness profile of galaxy cluster than the isothermal f3 model, except for outermost regions. 
Since the disagreement between our model and observed X-ray profiles in the outermost region 
is considered to be due to the employment of the King model, which describes a completely 
isolated system whose matter density becomes zero at a finete radius, it will disappear if the 
infall of matter from the outer (intercluster) space is taken into account. 



4. Discussion 

As shown above, our model is well fitted to the X-ray surface brightness profile except for 
outermost region. The fit is satisfactory even when the fit by the isothermal (3 model remains 
a large central excess (A 383 and A 2390). It should be noted here that our new model is well 
fitted to observed results for all of the sample clusters of galaxies (including both cD and non- 
cD types). The agreement between M pro j and Mi cns also looks like complete. The central mass 
concentration in our model is high enough to resolve the discrepancy between X-ray and lensing 
mass measurements. The ratio of the projected mass obtained using our model to the lensing 
mass is ranging from 0.92 to 1.06 for the three sample clusters, while the ratio obtained using 
the isothermal model ranges from 0.46 to 0.94 (see Table 2). This shows the predominance 
of our model over the isothermal (3 model. The agreement between M pro j and Mi cns would 
be improved if we take into account the infall of matter from the intercluster space. Since 
our model is based on firm physical foundations, these facts suggest that most of the 'central 
excess' of the X-ray surface brightness obtained so far would probably be no more than a ghost 
embossed by the use of inadequate model function. Therefore the central excess should not be 
directly linked to the cooling flow. 

Table 2 gives (3 = 0.53 ± 0.08 and 0.63 ± 0.02 for our model and the isothermal (3 model, 
respectively. If the velocity dispersion of intracluster dark matter is isotropic and if the in- 
tracluster gas has experienced no non-gravitational heating, (3 should not be less than unity 
because of the intracluster gas cooling by X-ray emissions (see Equation (6)). Therefore, the 



7 



values of f3 (our model) ~ 0.5 obtained here support the theory of non-gravitational intracluster 
gas heatings comparable to those reached through gravitational collapse (Kaiser 1966; Evrard 
& Henry 1991; Ponman, Cannon & Navarro 1999). The large value of (3 (our model) for A 2163 
is supposed to show the existence of some collective motions of intracluster gas. Indeed this 
cluster is not relaxed but just in the process of formation (Squires et al. 1997). 

So far we estimated the sizes and masses of sample clusters by assuming an Einstein-de 
Sitter (Qyi = 1.0, — 0.0) model of the universe for the convenience of the comparison of our 
results with other previous works. Here we consider the more realistic Qu = 0.3 and fl\ = 0.7 
cosmology too. Because, the r-dependence of the infalling matter density is sensitive to the 
value of Qm- Gunn & Gott III (1972) have calculated the infall of matter into the Coma cluster 
using the top hat model. Their calculations show that the infalling matter density varies as 
p oc r~ L67 for Q M = 2.0, oc r~ L55 for Q M = 0.5, and oc r~ L5 for Q M = 0.2. Accordingly 
we adopt p oc r -L52 for Qm = 0.3, which gives the integral of the square of infalling matter 
density along the line of sight varying as R l ~ 2xl - b2 — R~ 2m . It is noted here that this in- 
dependence will be valid for both of dark and baryonic matter particles, because the influence 
of gas pressure is negligible at the outermost region. In order to see the actual independence 
at large R (outer region) it is convenient to use the results of fit using the ispthermal f3 model, 
which gives S(R) oc R l ~ 6f3 there. The exponent -2.04 of R corresponds to the value of (3 = 0.51 
of the isothermal (3 model. This value is less than the observed mean value of f3 (isothermal) 
obtained here for the sample clusters and the mean value for 33 clusters (by the single (3 fit), 
< (3 >= 0.61 ±0.09 obtained by Xue & Wu (2000). These facts show that the matter density at 
cluster's outer region decreases with r more rapidly than estimated by the top hat model. This 
implies that the matter density in the space surrounding each cluster is fairly lower than the 
theoretically estimated intercluster mean matter density, probably due to the infall of matter 
towards each cluster over a long time. The recent discovery of a completely empty void of 
~ 140 Mpc radius (Rudnick et al. 2007) seems to support this point of view. 

As was mentioned in Subsection 2.1, Expression (9) is available even if A, T g andcr r are 
r-dependent so long as Expression (10) holds. In such a case, however, (3 would not have 
a constant value throughout a cluster, because T g and a r are expected to have different in- 
dependences in general. Now, from the Chandra or XMM-Newton observations, all of the 
sample clusters are found to have the gas temperature descent in the central region. The rates 
of temperature descent are respectively ~ 20% for A 383 (Furuzawa 2002), ~ 30% for A 2163 
(Markevitch & Vikhlinin 2001) and ~ 50% for A 2390 (Allen, Ettori & Fabian 2001). For each 
of these sample clusters, however, a very good fit of our model to the observed X-ray profile is 
obtained for r-independent values of T g and (3. Probably, the small change in the X-ray surface 
brightness caused by the gas temperature descent and the contributions of cD galaxy (when it 
exists) in the central region might be well compensated by a little but suitable adjustment for 
the best-fit values of an d (3. The parameter values compiled in Table 2 should therefore be 
understood as averaged over the whole system of each cluster. 

According to an attempt to derive the intracluster gas distribution embedded in a dark 
matter halo with an NFW density profile (Makino, Sasaki & Suto 1998), the obtained gas 
density profile is well approximated by the conventional (3 model (see also Suto, Sasaki & 
Makino 1998; Wu & Chiueh 2001; Wu & Xue 2000). In contrast to our model, however, some 
extra component is needed in case of the NFW dark halo to reproduce the central excess in 
the X-ray surface brightness distribution. The basic difference between our model and the 
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models based on the NFW model is focused on whether or not the central cusp exists for the 
dark matter distribution. This problem is still unsettled, and more detailed theoretical and 
observational studies are needed to bring this problem to an end. 

Aguirre, Schaye & Quataert (2001) and Sanders (2003) have tried to reproduce the observed 
distribution of intracluster gas on the basis of the modified Newtonian dynamics (MOND, Mil- 
grom 1983a, b, c) and showed that some additional dark matter component is still needed to get 
a good agreement between the modelled and observed results, especially in the central region of 
cluster. That is, the problem of the enhanced central mass concentration in clusters of galaxies 
is in principle separated from that concerning the validity of MOND. 

5. Conclusions 

A new model of intracluster gas distribution is established here. It is a new type of (3 model 
well fitted to the X-ray surface brightness profile except for outermost region of each cluster. 
The little gap between the modelled and observed X-ray surface profiles would disappear if the 
infall of matter from the intercluster space is taken into account. The fitting is satisfactory 
even when the X-ray surface brightness profile fitted to the isothermal (3 model remains a large 
central excess. The enhanced central mass concentration in our model is almost high enough to 
resolve the discrepancy between the X-ray and the gravitational lensing mass measurements. A 
little sink of the modelled X-ray surface brightness at the outermost region in our model would 
disappear if the infall of matter from the intercluster space is taken into account. These facts 
suggests that the central excess might be no more than a ghost embossed by an inadequate 
fitting and that the central excess should not be directly connected to the existence of cooling 
flow. 

The values of (3 (our model) ~ 0.5 obtained here support the theory of non-gravitational 
intracluster gas heatings comparable to those reached through gravitational collapse. On the 
other hand, the value of (3 ~ 0.6 obtained for the isothermal (3 model is explained as showing 
the scarcity of matter in the space surrounding each cluster, probably due to the infall of matter 
towards it over a long time. 
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Table 1. Redshifts, gas temperatures, arc radii, and 
lensing masses of sample clusters of galaxies. 



Cluster 


z 


T 


Rare 


Af lcns 






(keV) 


(kpc) 


(10 13 Af Q ) 


A 383 


0.188 


4.0 


65 




A2163 


0.201 


14 6 +0 - 9 


66 


2.4±0.8 


A 2390 


0.231 


ll.lil.O 


177 


16 ±2 



Table 2. Best-fit parameter values for our model and isothermal (3 model, and the resultant projected mass and 
projectcd-to-lcnsing mass ratio within the arc radius of each cluster. 



Cluster 




P 

(our model) 


P 

(isothermal) 


TO 

(kpc) 


a 

(kpc) 


M proj 
(1O 13 M ) 


MlT,proj 

(1O 13 M ) 


M proj 
M lcns 


MlT.proj 

M lens 


A 383 


s.o±8:S 


0.54 ±0.02 


0.65 ±0.03 


36 ±3 


93 ± 12 


3.8 ±0.2 


1.7 ±0.2 


0.95 ±0.31 


0.46 ±0.15 


A 2163 


7-4t°; 2 


0.60 ±0.01 


0.62 ±0.02 


310t*j 


305 ± 19 


2.2 ±0.2 


2.2 ±0.2 


0.92 ±0.34 


0.92 ±0.34 


A 2390 


8.2 ±0.2 


0.45 ±0.02 


0.62 ±0.01 


39+i° 


183 ±7 


17 ±2 


15 ±1 


1.06 ±0.17 


0.94 ±0.14 
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Figure 1: Simulated X-ray surface brightness plofiles for our model in case of (3 = 1.0: from 
top to bottom the central potentials of these models satisfy \l/o — 12, 9, 6, and 3. 
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Figure 2: X-ray surface brightness plofile of A 383 (Smith et al. 2001) fitted by our model (solid 
line) and the isothermal (3 model (dushed line). 
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Figure 3: X-ray surface brightness plofile of A 2163 (Elbaz, Arnaud & Bohringer 1995) fitted 
by our model (solid line) and the isothermal (3 model (dushed line). 
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Figure 4: X-ray surface brightness plofile of A 2390 (Bohringer et al. 1998) fitted by our model 
(solid line) and the isothermal (3 model (dushed line). 
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